Electronic transitions of single silicon vacancy centers in the near-infrared spectral 
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Photoluminescence (PL) spectra of single silicon vacancy (SiV) centers in diamond frequently 
feature very narrow room temperature PL lines in the near-infrared (NIR) spectral region, mostly 
between 820 nm and 840 nm, in addition to the well known zero-phonon-line (ZPL) at approx. 738 
nm [E. Neu et al., Phys. Rev. B 84, 205211 (2011)]. We here exemplarily prove for a single SiV 
center that this NIR PL is due to an additional purely electronic transition (ZPL). For the NIR 
line at 822.7 nm, we find a room temperature linewidth of 1.4 nm (2.6 meV). The line saturates 
at similar excitation power as the ZPL. ZPL and NIR line exhibit identical polarization properties. 
Cross-correlation measurements between the ZPL and the NIR line reveal anti-correlated emission 
and prove that the lines originate from a single SiV center, furthermore indicating a fast switching 
between the transitions (0.7 ns). g'^' auto-correlation measurements exclude that the NIR line is a 
vibronic sideband or that it arises due to a transition from/to a meta-stable (shelving) state. 
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INTRODUCTION 



In recent years, color centers in diamond have been in- 
vestigated as solid state light emitters especially in the 
context of single photon sources (e.g., Refs.[ll-ll) or, if in- 
corporated into nanodiamonds, as fluorescence labels for 
in vivo imaging (e.g., Refs. l6|-|8|). For these applications, 
color centers in diamond stand out owing to feasible room 
temperature operation as well as a high photostability. 
For applications in in vivo imaging, biocompatibility^ of 
the nanodiamonds and a large variety of possible surface 
functionalizations^ add to the advantages. The majority 
of previous experiments utilizes the well known nitrogen 
vacancy (NV) color center in diamond. Nevertheless, NV 
centers suffer from a major disadvantage, namely their 
broad room temperature emission bandwidth of about 
100 nm induced by strong electron-phonon coupling.^ 
Recently, silicon vacancy (SiV) centers emerged as a 
promising alternative j^!^!^ In contrast to NV centers, 
SiV centers enable narrow bandwidth room temperature 
single photon emission with a linewidth down to 0.7 nm 
and a high brightness (up to 5 Mcps).^ The narrow band- 
width results from low electron-phonon coupling leading 
to a significant concentration (exceeding 70%) of the lu- 
minescence in the narrow, purely electronic transition, 
i.e., the zero-phonon-line (ZPL), at approx. 738 nm. Fur- 
thermore, single SiV centers enable fully linearly polar- 
ized single photon emissionjiS which is advantageous for 
applications in quantum cryptographj*ii or in the fre- 
quency conversion of single photons SiV centers can 
be efficiently produced in situ by incorporation of Si 
impurities during chemical vapor deposition (CVD) of 
diamond. ^'^'^^ Furthermore, the production of colloidal 
solutions of fluorescent nanodiamonds containing SiV 
centers from polycrystalline CVD diamond films has been 
demonstrated. Thus, SiV centers are also promising can- 
didates for fluorescence labels especially due to feasible 
excitation with red laser light as well as their narrow 



emission in the red/near-infrarcd spectral rcgioui^ii^ For 
single photon emission as well as fluorescence labels, a 
further extension of color center emission into the longer 
wavelength near- infrared spectral range is of interest. De- 
spite the promising applications introduced above, sev- 
eral fundamental characteristics of the SiV complex are 
still unsettled: The charge state of the SiV center respon- 
sible for the 738 nm ZPL as well as the spatial symmetry 
of the complex are still under debateJ^iii In addition, 
a generally accepted suggestion for the electronic level 
scheme of the SiV center is still lacking. 

In this work, we aim at further elucidating the elec- 
tronic level scheme as well as the possible electronic tran- 
sitions of single SiV centers. In Ref. 03 we found that 
the photoluminescence (PL) spectra of single SiV cen- 
ters frequently feature very narrow room temperature PL 
lines in the near-infrared spectral region, mostly between 
820 nm and 840 nm. The present work is dedicated to 
an extensive investigation of these lines which we will 
term near-infrared (NIR) line(s). The NIR lines might 
arise from narrow vibronic sidebands or due to additional 
purely electronic transitions of the SiV centers. In previ- 
ous work, Sittas et alJ^ found a significant spectral nar- 
rowing upon cooling for three spectral features (776 nm, 
797 nm and 812 nm) in the sideband region of SiV cen- 
ter ensemble PL. Thus, they attributed these features 
to purely electronic rather than vibronic transitions. In 
earlier investigations, a three level model for single SiV 
centers has been established (see Refs. d and!3)- R in- 
cludes an excited and ground state giving rise to the ZPL 
transition and a third meta-stable level, termed shelving 
state. In the framework of this model, one might suspect 
that the NIR lines are due to a fluorescent transition from 
the excited state to the shelving state or due to a transi- 
tion from the shelving state to the ground state. In Ref. 
19, absorption measurements of single crystal CVD di- 
amond samples containing high densities of SiV centers 
are reported: The authors found absorption lines in a 
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wavelength range (830-860 nm) similar to the wavelength 
range of the NIR PL we observe here. The detection of 
a transition in absorption indicates that it involves the 
ground state of the system and only for purely electronic 
transitions, absorption and emission occur at the same 
wavelength. One thus might suspect that the NIR PL 
is attributed to an electronic transition from the shelv- 
ing state to the ground state. Similarly, for NV cen- 
ters, radiative transitions between shelving states (sin- 
glet states) have been reported at 1046 nm^S However, 
in the context of single color center spectroscopy, one has 
to make sure that the sharp PL lines do not belong to 
other color centers incorporated together with the SiV 
centers during CVD growth (e.g., nickel-related centers 
emitting between 806 nm and 820 nm).^^ In this work, 
we prove exemplarily for a single SiV center with pro- 
nounced NIR emission that the NIR line arises from an 
additional, independent purely electronic transition that 
cannot be explained within the established three level 
model. Section [n] introduces the investigated sample as 
well as the experimental setup. In Sec. lIIIi we discuss the 
spectrum of the SiV PL at room temperature. In Sec. lIVl 
we investigate the saturation curves of the ZPL and NIR 
PL. The polarization properties will be addressed in Sec. 
rvl In Sec. IVIl we present intensity auto-correlation g^^^^ 
measurements of the ZPL and the additional NIR line 
as well as cross-correlation measurements of both lines. 
The measurements prove single photon emission and re- 
veal the population dynamics of the SiV center. The 
cross-correlation measurements prove that the NIR line 
originates from the same color center as the well known 
ZPL. In Sec. lVIIl we discuss the measurements and iden- 
tify the NIR line as an independent electronic transition 
not connected to the shelving state. 



II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

The investigated SiV centers are contained in 
heteroepitaxial diamond films synthesized by mi- 
crowave plasma chemical vapor deposition (MPCVD) on 
Ir/YSZ/Si(001) substrates. The preparation of these 
substrates comprises the pulsed laser deposition (PLD) of 
an yttria-stabilized zirconia (YSZ) buffer layer followed 
by the e-beam evaporation of a single crystal Ir film on 
a 4° off-axis Si(OOl) wafer (for further details see Ref. 
[H). The mosaic spread of the Ir film is 0.15-0.2°. For 
the roughness we measure a value of <1 nm RMS. The 
bias enhanced nucleation (BEN) procedure is applied in 
a MPCVD reactor in order to generate heteroepitaxial 
diamond nuclei. During BEN the gas mixture in the re- 
actor consists of 3% CH4 in H2 at a pressure of 40 mbar. 
A plasma discharge is ignited by feeding in 2000 W of 
microwave power. A negative bias voltage of —300 V is 
applied to the Ir/YSZ/Si(001) substrate. After the BEN 
step, growth is performed without bias for 20 min at a 
substrate temperature of 800 °C and a reduced methane 



concentration of 0.5%. Scanning electron micrographs 
show a closed diamond layer. Its thickness is 90 nm as 
determined by ellipsometry measurements. SiV centers 
are created in situ with low density as a result of plasma 
etching of the Si substrates and subsequent incorporation 
of Si into the growing diamond)^ 

We investigate single SiV centers using confocal laser 
fluorescence microscopy. Optical excitation of the color 
centers is performed using a cw frequency-doubled, 
diode-pumped solid state laser providing an excitation 
wavelength of 671 nm corresponding to a photon energy 
of 1.85 eV. The excitation wavelength coincides with 
absorption bands of the SiV centeri^ Furthermore, it 
avoids direct photoionization of the SiV centers as the 
optical ionization threshold, i.e., the energetic distance 
of the ground state of the SiV center to the conduction 
band edge, is 2.05 eV (605 nm)i2^ The laser is focussed 
onto the sample by a high numerical aperture microscope 
objective (Olympus, magnification lOOx, NA 0.8). The 
fiuorescence is collected using the same objective and sep- 
arated from reflected laser light by a dichroic mirror and 
dielectric longpass fllters. For correlation measurements, 
we employ a Hanbury Brown Twiss (HBT) setup with 
two avalanche photodiodes (APDs, Perkin Elmer SPCM 
AQRH-14) featuring a typical quantum efficiency of ap- 
prox. 65% at 740 nm and 50% at 820 nm. Dielectric 
bandpass filters can be inserted in front of each APD to 
select the investigated luminescence lines (730-750 nm 
for the ZPL, 815-825 nm for the NIR line). Using this ex- 
perimental configuration enables the simultaneous mea- 
surement of saturation curves for the ZPL and the NIR 
line as well as intensity cross-correlation measurements. 
Photon arrival times are recorded with a fast timing elec- 
tronics (Pico Quant, Pico Harp, timing resolution of elec- 
tronics 4 ps, APD timing jitter 354 ps). We use this data 
to calculate the correlation functions. To investigate the 
polarized absorption of single SiV centers, the (linear) 
excitation polarization is rotated using a half-wave plate. 
The polarization of the emitted PL is investigated using 
a linear polarization analyzer. To investigate the spec- 
tra of the color center PL, we use a grating spectrom- 
eter (Horiba Jobin Yvon, iHr550). A grating with 600 
grooves/mm enables a resolution of approx. 0.18 nm. All 
experiments were performed at room temperature. 



III. SPECTROSCOPIC CHARACTERIZATION 

Figure [1] displays the spectrum of a single SiV cen- 
ter located in the (001) heteroepitaxial diamond film de- 
scribed in Secini The ZPL is clearly visible at 739.1 nm, 
a Lorentzian fit yields a linewidth of 0.9 nm (2.0 meV). 
The peak wavelength and the width of the ZPL is in ac- 
cordance with previous observations for single SiV cen- 
ters in randomly oriented nanodiamonds as well as het- 
eroepitaxial nanoislands.^-'^f We point out that we ob- 
serve a considerable inhomogeneous spread of the ZPL 
wavelengths of individual SiV centers in the investigated 
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heteroepitaxial diamond film in accordance with our pre- 
vious investigations:^*^ During the course of our investi- 
gations, single emitters with wavelengths ranging from 
735.3-739.2 nm have been observed. The inhomogeneous 
spread of ZPL peak positions is related to micro-stress 
fields at the locations of the individual color centers. 
Different stress sources are crucial in heteroepitaxial di- 
amond, e.g., stress due to thermal expansion mismatch 
of substrate and diamond as well as growth stress de- 
veloping between merging diamond grains. Since the re- 
lationship between the energetic shift of the SiV center 
ZPL and the local stress state is not known in detail, 
we cannot use the measured shifts to estimate the stress 
amplitude (for a detailed discussion see Ref. fiol). 

Figure [Ijb) shows the sideband region of the spec- 
trum in detail. As apparent from Fig. [Ijb), a single 
line dominates the spectrum in the wavelength range 
A > 800 nm. A Lorentzian fit to the narrow NIR line 
yields a peak position of 822.7 nm and a linewidth of 
1.4 nm (2.6 meV) [see inset of Fig. Ufa)]. Thus, the NIR 
line has a linewidth comparable to the ZPL. The energy 
difference between the ZPL and the NIR line is 170 meV, 
close to the energy of the Raman phonon mode in dia- 
mond (165 meV, see, e.g., Ref. ■26!). However, due to 
the different spectral widths as well as the results of the 
g^-^^ measurements discussed below, we exclude that the 
822.7 nm NIR line originates from Raman scattering of 
the 739.1 nm ZPL. Furthermore, the comparably high in- 
tensity of the NIR line excludes Raman scattered light. 
For single SiV centers (see Ref. 10), PL spectra often 
feature multiple NIR lines with wavelengths in the range 
from 820 nm to 840 nm, additionally excluding Raman 
scattered light. In low temperature experiments, we find 
a spectral narrowing of NIR lines down to 0.3 nm (0.56 
meV) at 5 K accompanied by a blue shift of 1.6 nm. In 
accordance with Ref. 0, this observation supports the 
identification of the NIR lines as electronic transitions. 
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FIG. 1. PL spectrum of a single SiV center. All spectra have 
been normalized to the ZPL. (a) Overview spectrum revealing 
a bright ZPL with low sideband contributions. The solid lines 
give Lorentzian fits to the data. The inset shows the dominant 
narrow PL line in the region A > 800 nm. (b) Zoomed image 
of the spectral region revealing the vibronic sidebands. 

The spectrum given in Fig. [T] does not straightfor- 
wardly allow for a comparison of the relative intensities 
of the emission lines as the detection efficiency varies over 



the observed wavelength range. Using the manufacturer 
supplied efficiency data for the spectrometer grating as 
well as the employed CCD detector, we find that the effi- 
ciency at 820 nm amounts to roughly 75% of the efficiency 
at 740 nm. Furthermore, the confocal setup has been op- 
timized for the detection of fluorescence at a wavelength 
of 740 nm, thus other wavelengths might not be optimally 
mapped. A more precise determination and comparison 
of the brightness of the emission lines is performed in Sec. 
IIVI using the single photon rates obtained with the APDs 
of the HBT setup. With this method, the setup can be 
individually optimized for the different emission lines. 



IV. SATURATION MEASUREMENTS 

As apparent from Fig. [Ija), the 822.7 nm line is sig- 
nificantly weaker than the ZPL at 739.1 nm. To quan- 
tify the PL intensity ratio, we perform a simultaneous 
measurement of the excitation power dependent photon 
count rates I{P) on these lines as described in Sec. |lTl 
Figure [2] displays the saturation curves measured on the 
ZPL and the 822.7 nm line. The fluorescence count rate 
I{P) obtained for a single color center is described by 



I{P) = /» 



P + Psat 



(1) 



if no linearly rising background luminescence is present. 
Psat is the saturation power, the maximum obtain- 
able photon rate. Fits using Eq. ([T|) describe the mea- 
sured data very well (see Fig. [5]), indicating a negligible 
contribution of background luminescence. For the ZPL, 
we find Psat = 61 ± 4 and = (55 ± 2) x 10^ cps. 
For the 822.7 nm line, we find Psat = 75 ± 4 and 
loo — (12.8 ± 0.4) X 10^ cps. The very similar saturation 
powers obtained for the two lines can be considered as 
first evidence that the PL is excited via the same absorp- 
tive transition. This in turn indicates that the lines stem 
from the same color center. Note that the count rate for 
the NIR line has been corrected for the lower detection 
efficiency of the APD at 820 nm as well as a lower band- 
pass filter transmission. Comparing the photon count 
rates, the ZPL intensity amounts to 4.3 times the inten- 
sity of the 822.7 nm line (the non-corrected count rate of 
the ZPL is about an order of magnitude higher than the 
NIR count rate). We point out that the NIR as well as 
ZPL emission were perfectly photostable: We did not de- 
tect any fluorescence intermittence during the course of 
measurements presented in this paper (checked for time 
intervals down to 1 ms) . The maximum excitation power 
used amounts to three times the saturation power Psat, 
proving photostability also for saturated excitation. 



V. POLARIZATION MEASUREMENTS 

Further information about luminescent transitions can 
be obtained using polarization spectroscopy.^'' First, we 
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FIG. 2. Saturation curves measured for the single SiV cen- 
ter. Filled squares (open circles) represent the datapoints, the 
solid line (dashed line) represents the fit using Eq. ([T]) for the 
ZPL (NIR line). The measurement for the NIR line has been 
corrected for the lower detection efficiency. 
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FIG. 3. Polarization measurements for a single SiV center: 
(a) measurements for the 822.7 nm line, (b) measurements 
for the 739.1 nm ZPL. Filled squares (open circles) represent 
the measured data of polarized absorption (emission polar- 
ization). Solid lines (dashed lines) give the sinusoidal fits to 
the absorption (emission) data. To enable comparison of the 
data, all measurements have been normalized to one. 



measure the polarized absorption of the SiV center by 
rotating the excitation polarization and simultaneously 
detecting the PL either on the 739.1 nm ZPL or on the 
822.7 nm line. Second, we measure the linear polariza- 
tion degree of the emitted light. 

Figure [3fa) [Fig. |31Ib)] displays the measurements for 
the 822.7 nm [739.1 nm] line in a polar plot (details see 
figure caption). The polarized absorption, i.e., the polar- 
ization dependence of the excitation efficiency, indicates 
the preferential absorption of linearly polarized light, 
consistent with a single absorption dipole, with a polar- 
ization direction of 89° for ZPL and NIR line. For light 
perpendicular to that direction, effectively no absorption 
takes place. We emphasize that polarized absorption only 
addresses the dipole component in the sample plane, i.e., 
the plane perpendicular to the excitation laser propaga- 
tion direction.— Thus, the two lines arise from an absorp- 
tive transition with the same in plane dipole moment. 
Note that an orientation of or 90° corresponds to an 
orientation along the (110) crystal directions in the (001) 
plane. Thus the observed orientation of the present in- 
dividual SiV center is in accordance with our previous 
observations in Ref. [l3 which were obtained using the 
ZPL emission. Together with the observation of a sim- 
ilar saturation power, this result supports the interpre- 
tation that the lines are excited via the same absorptive 
transition and thus arise from the same emitter. 

The PL shares a common linear polarization direction 
within the experimental error of approx. ±1° (84.0° for 
the 822.7 nm line and 85.6° for the 739.1 nm ZPL). Fig- 
ure[3l^b) shows a visibility close to 100%, while Fig.ISja) 
suggests a lower visibility. However, the reduced contrast 
is due to technical reasons, i.e., a reduced performance 
of the employed polarization analyzer in the NIR spec- 
tral range. Background fluorescence has been subtracted 
for both measurements; however, due to a spatially fluc- 
tuating background in the vicinity of the color center, 
this procedure may introduce an error in the visibility of 



5-10%. The observation of a common polarization direc- 
tion also indicates a common orientation of the emission 
dipoles in the sample plane. This result might be a hint 
that the lines originate from the same emitter. For the 
ZPL and the NIR line, polarized absorption and emis- 
sion yield an almost parallel orientation of the emission 
and absorption dipoles, respectively, (see Fig. [3]) in ac- 
cordance with previous findings for the ZPL of single SiV 
centers." 

The three-dimensional dipole orientation can only be 
retrieved using statistics for a large number of emit- 
ters that reveal the different equivalent directions (see, 
e.g., Ref. [lOl) . The investigation of the three-dimensional 
dipole orientation for the NIR transition is beyond the 
scope of this work. 



VI. INTENSITY CORRELATION 
MEASUREMENTS 

In this section, we investigate the luminescent transi- 
tions by measuring the intensity auto-correlation g^"^^ as 
well as the cross-correlation girliss of the PL lines, g*^^^ 
measurements can be used to prove the single photon 
nature of the PL as well as to analyze the population 
dynamics of the color center under investigation4i22ii^ 

(2) 

gcross measurements can prove that multiple fluorescence 
lines originate from the same emitter as demonstrated in 
Ref. 113 for the PL from different charge states of an NV 
center. Figure 0] summarizes the excitation power depen- 
dent intensity correlation measurements performed for 
the single SiV center. 

Figure[4l[a)-|-(b) displays the g^^-* measurements for the 

(2) 

739.1 nm ZPL (g^pi) revealing a distinct antibunching 
close to zero delay. For intermediate delays r, Qzpl 
exceeds one (bunching). To describe the intensity de- 
pendent g^pL functions, we use an extended three level 
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FIG. 4. Intensity correlation measurements (normalized as- 
suming g'^^ (r) — 1 for long delay times r). Note that adjacent 
g^^^ functions have been shifted for clarity [(a)-)~(c): adjacent 
jr'^-' functions shifted by 1, (b): adjacent g*-^' functions shifted 
by 0.5]: (a) gg^^ (739.1 nm line) 0.62 Psat - 2.92 Psat, (b) 
g'il^ 0.08 Psat - 0.37 Psat, (c) (822.7 nm line) with in- 

creasing excitation power (0.32 Psat, 0.81 Psat, 1.49 Psat), (d) 
intensity cross-correlation Qcross of 739.1 nm and 822.7 nm 
line (approx. 0.3 Psat)- (further explanation see text). Solid 
lines represent fits using Eq. ^ convoluted with the instru- 
ment response function of the HBT setup. 



scheme depicted in Fig. [T] The transition 2 — >■ 1 marks 
the ZPL. If the color center resides in state 3 (shelving 
state), no photons on the ZPL are emitted^ To obtain 
the (7^^-* function, one solves the rate equations for the 
populations N^: 



dNi 
"dT 
dN2 
dt 

<m 

dt 



= TVafei - Niki2 + A^afcai (2) 
= -iVafcai + Niki2 - N2k23 (3) 
= 7V2/C23 - Nsksi (4) 



We assume that the system is in the ground state at 
time zero and that the sum of the populations equals 
one. g^p^ is given by j^f^^^X'oo) (Ref. I2I) resulting in 



5zplW = !-(! + a) 



(5) 



The parameters a, ti and T2 are given byi^ 

ri,2= 2/{A±VA^~AB) (6) 

A= ki2 + k2i + k23 + ksi (7) 

B =ki2k23 + ki2k3i + k2lk3i + fc23fc31 (8) 

a = .'rj"'}. (9) 



The parameter ti governs the antibunching for short de- 
lays T, while the parameter T2 governs the bunching at in- 
termediate timescales. The parameter a determines how 
pronounced the bunching is. To describe the excitation 
power dependence of the paranieters a, ti and T2, we use 
the model introduced in Ref. |4|, where a linear depen- 
dence of the pumping rate coefficient ki2 of the excita- 
tion power P is assumed (ki2 = aP). Furthermore, the 
color center can be re-excited from the shelving state via 
a fourth state (see Fig. [7]) leading to a power dependent 
de-shelving rate coefficient k^i with 



«31 



d-P 

P + c 



k" 



(10) 



fe3l(T2-Ti) 



where k°i is the power independent de-shelving coeffi- 
cient and c and d describe the power dependence of ^31. 
In this model, all rate parameters can be determined us- 
ing limiting values of the parameters a, ri and T2 at high 
and low power as well as a fit of the power dependent 
curves of these parameters (for details see Ref. Q)- 

We use Eq. ([5|) convoluted with the instrument re- 
sponse function of our HBT setup to fit the measured 
data. The instrument response function of the setup is 
obtained independently via measuring the intensity auto- 
correlation function of attenuated ultrafast laser pulses 
from a mode-locked titanium-sapphire laser (Tsunami 
Spectra Physics, nominal pulse duration 100 fs). The 
instrument response is well described by a Gaussian func- 
tion with a full width at half maximum of 830 ps. The 
analytical result of the convolution of Eq. ([S]) with the 
instrument response is used as fitting function. The fit- 
ted functions are fully consistent with the measured data, 
thus proving single photon emission with negligible back- 
ground contributions. The fits provide the parameters a, 
Ti and T2 , summarized in Fig. [5] We deduce the param- 
eters describing the color center in the above specified 
model with intensity dependent de-shelving: ^23 = 5.8 
MHz, A:2i = 1.16 GHz, k°^ = 0.27 MHz, d = 18.96 MHz, 
cr = 8.4 MHz//iW and c = 66.6 fiW. Thus, the rate coeffi- 
cient for the ZPL transition (^21) is the largest coefhcient, 
while the rate coefficient populating the shelving state 
(/c23) is significantly smaller. The de-shelving rate coeffi- 
cient at low power (k^^) is nearly an order of magnitude 
smaller than its high power limit (d + k^i). Figure [5] also 
displays the power dependent curves for ti, T2 and a ob- 
tained from this model as solid lines and the curves from 
a model with constant k^i as dashed lines (for details see 
Ref. 4). As apparent from Fig. [5l the increase of T2 at 
low excitation powers is overestimated by the intensity 
dependent de-shelving model, whereas it is not at all de- 
scribed using the model with constant de-shelving. The 
power dependence of ri and a is reasonably described in 
the model including intensity dependent de-shelving, fit- 
ting the data more closely than the model with constant 
k^i. These observations thus support the assumption of 
an intensity dependent de-shelving process. 

Figure Wic) displays the g^^^ auto-correlation mea- 

(2) 

surements of the 822.7 nm fluorescence (ffJv/fl)- Fewer 
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measurements have been performed as the measurement 
times are significantly increased due to the low count 
rates. The fits of gjv/i? El'^) show that the devi- 

ation from 5(^^(0) = is fully explained by the instru- 
ment response of our HBT setup which has been taken 
into account for the fits. As apparent from Figs. \M,c) 
and[5l g^^\ji displays only a very weak bunching. Fitting 
these g^^'^ functions analogously to the ones of the ZPL, 
the antibunching time constants Ti^^ are about a factor 

of 3 shorter than for g^^pj^ (see Fig. [SJ low power limit 
^NiR.M ^ Q 26 ns, = 0.86 ns). 



3600 
2400 
1200 


■ 


■ 

r-i 


. 739.1 nmg'" - 
822.7 nmg'" ; 
X Cross correlation 

■ 







1.00 
0.75 
0.50 


■ 

^^^^ 


■ 


^ ■ 


0.25 











0.48 






■ 


0.32 




■ 


■ 


0.16 
0.00 











50 100 150 200 

Excitation Power (|jW) 



FIG. 5. Parameters a, ri and T2 obtained for the g'^' mea- 
surements. Black solid lines give the intensity dependent 
curves for the ZPL g^^^ function in the intensity dependent 
de-shelving model. Gray, dashed curves are plotted using the 
model with constant de-shelving (explanation see text). 

In addition to the g^^-* auto-correlation measurements, 
we perform cl Qcr^oss 

cross-correlation measurement of the 

ZPL and NIR emission lines. If these lines originate from 

(2) 

the same emitting color center, a gcross measurement re- 
veals an antibunching: the emission is anti-correlated. If 
the lines originate from different emitters, one expects 
no correlation between the fiuorescence photons 1^ We 
point out that the gcross measurement does not enable 
to discriminate between lines from different charge states 
of the emitter or lines from different transitions (purely 

electronic or vibronic) of the same charge state. 
(2) 

gcross measurements have been used to prove the 
charge state conversion of a single NV center from NV^ 
to NV°i^ We here assume that the 822.7 nm line is not 
due to a second charge state, as SiV° emission has been 
reported at 946 nm.-^- Nevertheless, also the improbable 
situation of a single electron alternatively charging an 
SiV center (emitting at 739.1 nm) and another color cen- 



ter (emitting at 822.7 nm) would lead to anti-correlated 
emission. However, one might also expect to observe 
fiuorescence blinking in this case, which we do not find 
in our experiments. Additionally, the similar saturation 
curves and polarization measurements support the as- 
sumption that the lines arise from a single emitter. Fig- 

(2) (2) 

ure m^d) displays the measurement of gcross- gcross has 
been obtained using two methods. First, we calculate the 
cross-correlation by correlating each time tag (photon ar- 
rival time) recorded on one avalanche photodiode with all 
other time tags for the second diode [filled squares in Fig. 
lU^d)] analogously to the auto-correlation measurements. 
Second, we simulate a start-stop measurement by using 
each detection event on the NIR line as a start event 
and search for the next emission of a ZPL photon which 
serves as the stop event and determines the delay r [light 
gray line in Fig.Hl^d)]. The second method has previously 
been used to obtain cross-correlation functions for quan- 
tum dot emission lines. ■^^'"^^ We find strong antibunching, 
the curves calculated using the two methods are equal 
within experimental errors. The absence of an asymme- 
try in the start-stop measurement proves that there is 
no preferred time-ordering of the emission events of both 

lines (see also discussion Sec. IVII[) . The time constant 

(2) 

^cross ^0.7 ns is close to the values obtained for g^pj^ 

nv (2) 

[fit: black line in Fig.UJJd)]. In summary, the gcross mea- 
surement identifies the 822.7 nm PL as anti-correlated to 
the SiV ZPL. Thus, we conclude that the lines originate 
from the same emitter. 



VII. DISCUSSION OF THE CORRELATION 
MEASUREMENTS 

Using the g^"^^ measurements introduced above, we first 
deduce the origin of the 822.7 nm luminescence: First, 
we emphasize that the different time constants ri and 
^NiR exclude that the NIR line is a vibronic sideband: 
Vibronic sidebands originate from the same excited state 
as the ZPL; however, they end in vibrationally excited 
states of the ground state.-^^ These states decay within 
picoseconds to the vibrational ground state which is also 
the final state of the ZPL transitioui^ This very fast 
relaxation following the emission of a photon does not 
influence the measured g^^^ and thus the g^^^ function 
measured for sideband fluorescence is equal to the 17 
function measured for the ZPL fluorescence itself. 

Second, on the basis of the level scheme depicted in 
Fig. [7] (Level 1-4) used to model the ZPL emission dy- 
namics one might suspect that the 822.7 nm PL arises 
from the de-shelving transition 3 — > 1. A participation 
of the ground state is motivated by the observation of 
absorbing transitions at similar wavelengths in diamonds 
containing SiV ensembles. For such a situation, we cal- 
culate the g^^-' function for the light emitted on the tran- 



7 



(2) 

sition 3^-1 analogously to g^pL 



(2) 



Nsir oo) 



(11) 



as also here each photon detection event projects the sys- 
tem into the ground state and we can use the same initial 
condition as for the transition 2—^1. From the discus- 
sion of the ZPL emission dynamics, fcsi was found to be 
intensity dependent [Eq. dTU])]. For our model, we as- 
sume that the transition 3 — > 1 is radiative. In contrast, 
the intensity dependent de-shelving process can be in- 
terpreted as a non-radiative transition 3 — > 1 via state 
4. Both processes modify the population of state 3 and 

(2) 

thus g^J^i- Nevertheless, for each excitation power, the 
probability for a photon emission on the transition 3 — > 1 
after a delay r is proportional to the population N^It) 
and thus Eq. (ITT|) indeed describes the g^^'^ function. 

We now solve the coupled rate equations [Eqs. ©-(H))] 
using k2i, fc23, fcg]^, d, c and a obtained from the measure- 
ments of g^zPL"^ From the solution of the rate equations, 
we calculate the normalized occupation for states 3 and 
2 for different excitation powers P. Figure [6] summarizes 

f 21 

the resulting normalized occupation of state 2 {gzpL) 

f2) 

well as of state 3 ((734.1), details see figure caption. Com- 
paring Figs. |4Ua)-|-(b) and[6l it is apparent that the sim- 
ulated g^"^^ functions reasonably agree with the measured 
9 ZPL functions, thus the parameters deduced using the 
intensity dependent de-shelving model appropriately de- 
scribe the observed SiV center. The simulated functions 
do not include a correction for the instrument response, 
thus they display g'^'^\Q) = 0, whereas the measured g*^^-' 

functions display a non- vanishing value for ^'■^•'(O). 

f 21 

As evident from Fig. [51 g^_^i also displays an anti- 
bunching; however, the corresponding time constant ^ 
is significantly longer than ri or t(^^^: Below saturation 

[Fig. [6ja) and (b)], g^J^i raises to a constant value of 

one within several 100 ns as the rate coefficients ^23, fcai 

(2) 

are much smaller than fc2i. Thus, the simulated 534.1 
is in stark contrast to gj^jp measured on the 822.7 nm 
line [see Fig. UJ^c)]. The (/^^■' measurements thus exclude 
that the 822.7 nm line arises due to a transition between 

the shelving state (state 3) and the ground state (state 
f 21 

1). Furthermore, the gj^jp measurements exclude that 
the NIR line arises due to the transition 2 — >■ 3: For a 
process involving state 2 as an intermediate state, t{^^^ 
should be at least as long as ti. 

We emphasize that the ZPL and NIR line cannot arise 
from a cascaded emission as the cross-correlation mea- 
surements do not reveal a bunching, thus we can exclude 
that one of the transitions populates the excited state 
of the other transitionJi Furthermore, the photon en- 
ergy of the excitation laser is 1.78 eV which is well below 
the sum of the energies for both transitions (1.68 eV and 
1.51 eV). Additionally, the linear dependence of the fluo- 
rescence on excitation power for both transitions excludes 
multi-photon processes for the excitation. 
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FIG. 6. Normalized occupation of state 3 and state 2, simu- 
lated using fc23 = 5.8 MHz, k2i = 1.16 GHz, = 0.27 MHz, 
d = 18.96 MHz and a = 8.4 MRz/fiW and c = 66.6 ixW. (a) 
P = 7 MW,(b) P = 35 fiW, (c) P = 140 /iW, (d) P = 210 
fiW. The excitation powers roughly correspond to O.lPsat, 
0.5 Psat, '2 Psat and SPsat- Insets show zoomed images. 
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NIR 
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FIG. 7. Schematic representation of the extended model in- 
cluding the NIR transition. The level scheme included in the 
dashed box indicates the level scheme for the intensity depen- 
dent de-shelving model. Length of the NIR transition and of 
the ZPL transition drawn to scale. 



(2) 

As the gcross measurements prove that the NIR line 
and the ZPL originate from the same emitter, we propose 
an extended level scheme including the NIR transition as 
depicted in Fig.[71 On the basis of the evidence presented 
above, we assume that the NIR line originates from an 
independent electronic transition, i.e., a second ZPL. We 
introduce the NIR transition 5 — ?> 1 as a third relaxation 
path into the level scheme (Fig. [7]). As g^^PL matches 
the g^"^^ function of a three level system, state 5 should 
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not accumulate population (as it is the case for state 3): 
This would induce deviations from the three level g*^^^ 
function. The short antibunching time constant t^^^' 
supports the assumption of a fast decay path. Thus, 
this pathway, similar to fast non-radiative decay paths, 
might not lead to a deviation from the g*^^^ function of 
a three level system. We assume that the relaxation to 
state 5 occurs from the intermediate pumping level in 
the off-resonant pumping process: For a process involv- 
ing state 2 as an intermediate state, t^^^ should be at 
least as long as ti , as the filling of state 5 from state 2 
would limit the antibunching time constant. The relax- 
ation rate coefficients from the pumping level are not di- 
rectly accessible in our measurements. We assume these 
rate coefficients significantly exceed any other coefficients 
in the system ensuring that the pumping level does not 
accumulate population. This assumption is justified, as 
the vibrationally excited states that are most probably 
responsible for the pumping transition used here feature 
very short lifetimes in the picosecond range The ratio 
of the relaxation rates from the pumping state to state 
5 and to state 2 should infiuence the branching ratio for 
the 740 nm ZPL and the NIR line. However, taking into 
account also a potentially differing quantum efficiency, 
we cannot use the intensity ratio of the lines to deduce 
the branching ratio. 

For a two level system at low excitation powers, ti in- 
dicates the lifetime of the excited state of the radiating 
transition assuming no longer lived states are populated 
in the excitation process. Thus, from the g^^^ measure- 
ments, one might expect a very bright emission from the 
822.7 nm line due to a short lifetime. However, the frac- 
tion of the population available for the 822.7 nm line is 
unsettled taking into account the unknown branching ra- 
tio of the transition from the pumping level to states 
2 and 5. In principle, the lower brightness might be 
linked to a more efficient quenching mechanism for the 
NIR luminescence resulting in a low quantum efficiency.'^*' 
Non-radiative quenching of low energy transitions has 
been reported for NV~ centers: The infrared emission 
at 1046 nm is four orders of magnitude weaker than the 
visible 637 nm emission despite the fact that 30% of the 
population is available for the infrared transitioni^ 

(2) 

In the following, we further analyze the Qcross measure- 
ments and infer the switching dynamics of ZPL and NIR 
line. Qcross reveals a fully symmetric antibunching. An 
asymmetric antibunching would indicate that the time 
elapsing between the emission of an NIR photon and a 
ZPL photon is different from the time elapsing for the re- 
verse order of events i^^'^^ In our model, the two lines are 
populated via the same pumping transition and share the 
same ground state, consistent with the observation of a 
symmetric cross-correlation. For the interpretation of the 
cross-correlation time constant rf"*'*, we consider Ref. 

IstI where f/irLs measurements between disjunct spectral 
windows of the PL transition of a single quantum dot 
have been introduced. The emission line is broadened 
by spectral diffusion, consequently a change of the emis- 



sion wavelength from one spectral window to the other 

(2) 

requires a spectral jump of the emitter, gcross exhibits 
an antibunching, where rf"''^ reveals the characteristic 
time Td for the spectral jumps. Interpreting our mea- 
surements analogously, rf""^ reveals the characteristic 
time for the SiV center to change its emission between 
the NIR and ZPL transitions. As rf"''* is comparable 
to the lifetime of the excited state (state 2), the changes 
occur in between successive excitation and emission cy- 
cles of the ZPL. This short characteristic time also sup- 
ports the assumption of an alternative decay path con- 
necting the pumping levels and the ground state of the 
center without involving the long lived shelving state: 
The longer time windows, in which the emitter resides in 
the shelving state (state 3) do not influence the measured 

(2) 

gcross function between the ZPL and the NIR lines if the 
level scheme in Fig. [7] is applicable. Additionally, gcross 
indicates that the SiV center changes between the two 
emission lines frequently instead of undergoing a larger 
number of transitions on cither the ZPL or the NIR line 
before changing its emission wavelength again. 



VIII. CONCLUSION 



Photoluminescence (PL) spectra of single SiV centers 
frequently feature very narrow room temperature PL 
lines in the near-infrared spectral region, mostly between 
820 nm and 840 nm.W In this work, we clearly prove for 
a single SiV center, that this PL (at 822.7 nm) is due 
to an additional electronic transition (besides the well 
known ZPL here at 739.1 nm). We find a linewidth of 
1.4 nm (2.6 meV) for the NIR line. The NIR line satu- 
rates at similar excitation power as the ZPL. However, 
despite a shorter excited state lifetime deduced from the 
g^^) measurements, it delivers a factor of four lower fiu- 

(2) 

orescence intensity. A gcross measurement between the 
ZPL and the NIR line reveals anti-correlated emission, 
thus proving that the emission originates from the same 
emitter. The short antibunching time constant in the 

(2) 

gcross measurements indicates a fast switching between 
the transitions, g^^^ auto-correlation measurements ex- 
clude that the NIR emission arises due to a transition 
from/to the shelving state of the SiV center. They fur- 
ther exclude that the NIR line is a vibronic sideband. 
Polarization measurements reveal preferential absorption 
of linearly polarized light for both transitions. The max- 
imum absorption is observed for the same polarization 
direction. The emitted fluorescence light shares a linear 
polarization direction. We interpret the NIR transition 
as an independent electronic transition populated from 
the same pumping levels as the ZPL transition. 



9 



ACKNOWLEDGMENTS 



We acknowledge funding by the DFG and the BMBF 
(EPHQUAM 01BL0903). 



* matthias.schreck@physik.uni-augsburg.dej 

t £hristoph. becher@physik.uni-saarland.de 

^ R. Brouri, A. Beveratos, J. Poizat, and P. Grangier, Opt. 
Lett. 25, 1294 (2000). 

^ C. Kurtsiefer, S. Mayer, P. Zarda, and H. Weinfurter, 
Phys. Rev. Lett. 85, 290 (2000) 

^ C. Wang, C. Kurtsiefer, H. Weinfurter, and B. Burchard, 
J. Phys. B: At. Mol. Opt. Phys. 39, 37 (2006). 

^ E. Neu, D. Steinmetz, J. Riedrich-Moller, S. Gsell, M. Fis- 
cher, M. Schreck, and C. Becher, New J. Phys. 13, 025012 
(2011). 

^ T. Gaebel, L Popa, A. Gruber, M. Domham, F. Jelezko, 
and J. Wachtrup, New J. Phys. 6, 98 (2004) 

^ Y.-R. Chang, H.-Y. Lee, K. Chen, C.-C. Chang, D.- 
S. Tsai, C.-C. Fu, T.-S. Lim, Y.-K. Tzeng, C- 
Y. Fan g, C.-C. Han, H.-C. Chang, and W. Fann, 
[Nature Nanotech. 3, 284 (2008) 

N. Mohan, C.-S. Chen, H.-H. Hsieh, Y.-C. Wu, and H.-C. 

Chang, Nano Lett. 10, 3692 (2010) 
® F. Neugart, A. Zappe, F. Jelez koVC. Tietz, J . P. Boudou, 

A. Krueger, and J. Wrachtrup, ,Nano Lett. 7, 3588 (2007) 
^ E. Neu, C. Arend, E. Gross, F. Guldner, C. Hepp, 

D. Steinmetz, E. Zscherpel, S. Ghodbane, H. Stern- 

schulte, D. Steinmiiller-Nethl, Y. Liang, A. Krueger, and 

C. Becher, Appl. Phys. Lett. 98, 243107 (2011) 
^° E. Neu, M. Fischer, S. Gseh, M . Schreck, and C. Becher, 

[Phys. Rev. B 84, 2052H (20il ) 

C. Bennett, G. Brassard, et al, Proceedings of IEEE In- 
ternational Conference on Computers, Systems and Signal 
Processing, Bangalore, India 175 (1984). 
S. Zaske, A. Lenhard, and C. Becher, Opt. Express 19, 
12825 (2011). 

J. Barjon, E. Rzepka, F. Jomard, J. Laroche, D. Ballutaud, 
T. Kociniewski, and J. Chevallier, Phys. Status Solidi A 
202, 2177 (2005). 

G. DoUinger, A. Bergmaier, C. Frey, M. Roesler, and 

H. Verhoeven, Diam. Relat. Mater. 4, 591 (1995). 

R. Weissleder and V. Ntziachristos, Nat. Med. 9, 123 
(2003). 

S. Moliver, Tech. Phys. 48, 1449 (2003). 

J. P. Gos s, R. Jones, S. J. Breuer, P. R. Briddon, and 

S. 6berg,'Phys. Rev. Lett. 77, 3041 (1996) 

G. Sittas, I. Kiflawi, H. Kanda, and P. Spear, Diam. Relat. 

Mater. 5, 866 (1996). 



U. F. S. D'Haenens- Johansson, A. M. Edmonds, B. L. 
Green, M. E. Newton, G. Davies, P. M. Martineau, 
R. U. A. Khan, and D. J. Twitchen, Phys. Rev. B 84, 
245208 (2011). 

L. J. Rogers, S. A rmstrong, M . J. Sellars, and N. B. Man- 
son, New J. Phys. 10, 10"302 4 (2008) 

I. Aharonovich, C. Zhou, A. Stacey, 
F. Treussart, J.-F. Roch, and S. Prawer, 

Appl. Phys. Lett. 93, 243112 (2008)[ 

E. Wu, V. Jacques, H. Zeng, P. Grangier, F. Treussart, 
and J.-R. Roch, Opt. Express 14, 1296 (2006). 
S. Gsell, M. Fischer, M. Schreck, and B. Stritzker, 
J. Cryst. Growth 311, 3731 (2009) 

K. lakoubovskii, G. Adriaenssens, N. Dogadkin, and 

A. Shiryaev, Diam. Relat. Mater. 10, 18 (2001). 

K. lakoubovskii and G. J. Adriaenssens, 



Phys. Rev. B 61, 10174 (2000) 



A. 



K. 



Ramdas, 



~K'. Solin and 
IPhys. Rev. B 1, 1687 (1970) 
A. Kaplyanskii, Opt. Spec 16, 329 (1963). 
T. Ha, T. Laurence, D. Chemla, and S. Weiss, J. Phys. 
Chem. B 103, 6839 (1999). 

I. Aharonovich, S. Castelletto, D. A. Simp- 
son, A. D. Greentree, and S. Prawer, 
Phys. Rev. A 81, 043813 (2010) 
^° T. Gaebel, M. Domhan, C. Wittmann, I. Popa, 

F. Jelezko, J. Rabeau, A. Greentree, S. Prawer, 
E. Trajkov, P. Hemmer, and J. Wrachtrup, 
lAppl. Phys. B: Lasers Opt. 82, 243 (2006)| 

A. Khaz, S. Faith, C. Becher, B. Gayral, W. V. Schoen- 
feld, P. M. Petroff, L. Zhang, E. Hu, and A. Imamoglu, 
jPhys. Rev. B 65, 16130 3 (2002) 

D. V. Regelman, U. Mizrahi, D. Gershoni, E. Ehren- 
freund, W. V. Schoenfeld, and P. M. Petroff, 
"Phys. Rev. Lett. 87, 257401 (2001) 
J. Walker, Rep. Prog. Phys. 42, 1605 (1979). 

B. Lounis and M. Orrit, Rep. Prog. Phys. 68, 1129 (2005). 
Solving of the differential equations (analytically) and cal- 
culating the limiting values 7V2,3(r — >■ oo) performed using 

Maplesoft ' s Maple 15. 

L. Rogers, [Physics Proced ia 3, 1557 (2010) | 10th Interna- 
tional Meeting On Hole Burning, Single Molecule, And 
Related Spectroscopies, 2009, Palm Cove, Australia. 

G. Sallen, A. Tribu, T. Aichele, R. Andre, L. Besombes, 

C. Bougerol, M. Richard, S. Tatarenko, K. Kheng, and 
J. Poizat, Nature Photon. 4, 696 (2010). 



33 
34 



